Scaling and interaction-assisted transport in graphene with one-dimensional defects 
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We analyze the scattering from one-dimensional defects in intrinsic graphene. The Coulomb 
repulsion between electrons is found to be able to induce singularities of such scattering at zero 
temperature as in one-dimensional conductors. In striking contrast to electrons in one space dimen- 
sion, however, repulsive interactions here can enhance transport. We present explicit calculations 
for the scattering from vector potentials that appear when strips of the material are under strain. 
There the predicted effects are exponentially large for strong scatterers. 

PACS numbers: 72.80.Vp, 71.45.Gm, 72.10.Fk, 73.63.-b 
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The Coulomb repulsion between electrons can have 
profound consequences for the scattering of electrons in 
conductors. In one-dimensional (ID) conductors it dra- 
matically suppresses the conductance through impurities, 
with a singularity at temperature T = 0. This is one of 
the hallmarks of the Luttinger liquid state of interacting 
electrons in ID |l|. A particularly intuitive understand- 
ing of this effect is due to Matveev, Yue, and Glazman 
[1, 01]. In their approach interactions cause extra elec- 
tron scattering from Friedel oscillations, which inhibits 
transport. Friedel oscillations are density modulations 
that form when electron waves incident on an impurity 
interfere with backscattered waves. In ID at T = the 
Friedel oscillation around an impurity is not integrable 
over space and the induced scattering amplitude diverges 
logarithmically in Born approximation. 

Friedel oscillations in two dimensions (2D) at nonzero 
chemical potential can be integrated and have much 
weaker effects for both, point-like and ID de- 

fects @, HI- Interactions play a more important role 
in graphene [il-HH at zero chemical potential ("intrin- 
sic graphene"), which forms a so-called marginal Fermi 
liquid 12j. There is, for instance, experimental evidence 
13] of a singular renormalization of the Fermi velocity 
12j,ll4|. Nevertheless, also in intrinsic graphene scatter- 
ing from point-like impurities of both, scalar and vector 
potential character, does not receive any singula r correc- 
tions from the Coulomb interaction 



In this Letter we show that the Coulomb interaction 
in intrinsic graphene can have a singular impact on scat- 
terers that are extended in one space dimension. A qual- 
itative analysis of the Friedel oscillations at such scat- 
terers reveals the reason: the Friedel oscillation due to 
an electron with wavevector k' scattering with ampli- 
tude r from a ID defect fills a strip of width ~ 1/k'. 
The corresponding exchange potential has support in the 
same strip, implying a factor ~ r/k' in the induced Born 
scattering amplitude Sr cx . A second factor 1/k' comes 
from the Coulomb interaction. For energy-independent r 
the exchange interaction with electrons in all filled states 
thus results in <5r cx oc r / fillcd cPk'/k' 2 . Since in intrinsic 
graphene at T — the entire valence band is filled, (5r ex 



is logarithmically divergent. This suggests that when the 
Coulomb interaction does not cause a Hartree potential it 
can have the same drastic effects for ID scatterers in 2D 
intrinsic graphene as for point defects in ID conductors. 

In contrast to the familiar situation in ID [l[, how- 
ever, we demonstrate that the Coulomb interaction in 
graphene can enhance electron transport. This intriguing 
effect is due to bound states that straddle defects. The 
exchange interaction of transport electrons with electrons 
in such bound states opens an additional transport path 
across the defect, which increases the electric current. 




FIG. 1: A vector potential A y in a narrow strip (dark blue) 
forms a scattering barrier in a sheet of otherwise ballistic 
graphene (light blue with contacts in green). A y may be in- 
duced by strain or, alternatively, by a pair of current-carrying 
wires (brown). 



We exemplify the above effects with the scattering bar- 
rier shown in Fig. [TJ formed by a vector potential in a 
narrow strip of an otherwise ballistic sheet of intrinsic 
graphene. Such vector potentials are induced by strain 
17, [l8| or by nearby electric currents, see Fig. [T] By 
symmetry vector potentials in graphene do not produce a 
Hartree potential. As suggested by the above argument, 
the Coulomb interaction therefore induces extra scatter- 
ing whose amplitude diverges at low T, when the elec- 
tron wavelength is much longer than the width a of the 
barrier and scattering is energy- independent. The effect 
turns out to be exponentially enhanced for strong barri- 
ers, where it strongly impacts transport even at moder- 
ate T. In addition, the scatterer of Fig. [1] hosts bound 
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states that enhance transport through electron exchange 
as motivated above. As a consequence, the barrier of 
Fig. [T] becomes entirely transparent at T = 0, in striking 
contrast to impurities in interacting ID conductors that 
completely block the electric current at T = 0. 
Model: In the setup of Fig. Q] a vector potential 
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[@(x) is the step function] is applied to a sheet of clean, 
intrinsic graphene that we first assume to be infinite, L — ► 
oo. The noninteracting Hamiltonian in valley 7 thus is 



iJ 7 = vcr-y ■ (p — eA) 



(2) 



Here, er 7 — (a x ,ja y ) is a vector of Pauli matrices, p is 
the electron momentum, v the electron velocity, and — e 
the electron charge. The above model has parity sym- 
metry, VHP = H, where Pip(x,y) = a y ^p(—x,y), and 
particle-hole (PH) symmetry, o~ z Ha z — —H. In intrinsic 
graphene both symmetries are also respected by electron- 
electron interactions [24| . Without restriction we assume 
X > 0. Transport through the barrier at x < is ob- 
tained through reflection on the x-axis. The above model 
has been analyzed in Ref. [l7| on a noninteracting level, 
where A suppresses linear transport exponentially in x if 
kT -C X v / a ( we se t ^ — !)• Below we study the effects of 
the Coulomb repulsion between electrons on this trans- 
port problem for short barriers (Aa ^> 1, but InAa > 1, 
where A is the bandwidth of graphene) . 

We assume that the graphene sheet is suspended or 
lying on an insulating substrate, leaving the Coulomb in- 
teraction unscreened with interaction potential Vc (q) = 
2irr s v/\q\ at interaction parameter f s and wavevector q 



[25j |. In the weakly interacting limit f s C 1 that we take 
henceforth the main effect of electron-electron interac- 
tions is additional elastic scattering, as in the thoroughly 
studied ID problem 0, El El- Also like in ID, such scat- 
tering is induced by the Hartree-Fock potential V HF due 
to Friedel oscillations at the barrier. In addition, A here 
induces bound states with an extra contribution to y HF . 

The potential V^ F is strongly constrained by the sym- 
metries of the setup. Most importantly, a Hartree po- 
tential is precluded by PH-symmetry. Also, the ex- 
change potential V cx is spin- and valley-diagonal. Trans- 
port in both valleys is identical since they are related 
as = U$HiU v , where U v = o~ x commutes with the 
current er x through the barrier [3]. We thus set 7 = 1. 

The compound scatterer composed of A and V cx is 
conveniently described by the transfer matrix of elec- 
trons at zero energy M — T(oo, — oo)| e= o Q. The V- 
and PH-symmetries of our model imply o~ y Mo y = M 
and o~ z Mcr z — M, respectively. Together with current 
conservation er x AfV x = M~ l they constrain M to take 
the form M = exp(xc z ). Low energy transport in the 
setup of Fig. [1] with weak many-body interactions is thus 
characterized by a single parameter x- This reduces our 
analysis to a calculation of x- At f s = we have x — X- 

First order in f s : To begin with, we calculate the in- 
teraction correction to x a t &rst order in f s . Computing 
the correction to M due to scattering from the exchange 
potential created by the noninteracting electron states in 
Born approximation we find in the limit a <C v/kT 



X = X~ r s s'mhxF(x) 
with the positive function 



In- 



kTa 



X + 0(l) (3) 



, , coth v(cosh3Y — 3 cosh y) arcsin(tanh x) + 2 cosh 2 y . 
F( X ) = ^ - cosh X smh 2 x- (4) 



The first order result Eq. (j3|) manifests itself in 
a non-analytic temperature-dependence of the conduc- 
tance through the barrier, which takes the form 

4e 2 kTWln2 / , , v\ , r , 

G = — 1 + smhxtanhx hitanh — J . (5) 

h 7tv V 2 / 

In G, the logarithmic temperature dependence of the in- 
teraction correction to x competes with the factor T due 
to the linear density of states. The predicted non-analytic 
scale-dependence of x is thus more easily observed in the 
normalized conductance G T = G/(G|^=o) [Ill- 
Due to its logarithmic T-dependence the first order 
correction to x may grow large even at f s <C 1- At the 
corresponding low temperatures, when f s \n(v/kTa) > 1, 



the above perturbative calculation looses validity, as in 
the ID case The origin of the logarithmic diver- 

gence here is, however, different. In ID the divergence 
stems from a Friedel oscillation density oc 1/x at T = 0. 
In the present 2D case the T — Friedel oscillation den- 
sity is cx 1/x 2 and integrable. But the non- locality of V ex 
gives rise to the same logarithmic divergencies as in ID. 

RG- analysis: We extend our analysis into the regime 
f s <C 1, but f s \n(v /kTa) > 1, by a renormalization group 
(RG) calculation that re-sums the perturbation series in 
the "leading logarithm" approximation 0, Q . To this end 
we successively integrate out shells of wavevectors [k 1 , bk'] 
with b > I in the order of decreasing k 1 <C 1/a. For 
each shell we compute the exchange potential V£,\ due to 
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the corresponding states. Thereafter we renormalize the 
transfer matrix for low energy electrons at wavevector k 
by V£,\. The resulting x characterizes the wavefunctions 
in the subsequent shell [k, bk]. The logarithmic scale- 
dependence of Sx, Eq. ([3]), allows us to assume k <C k' in 
this process. Introducing I = — In k'a we thus find that 



dl 



-r s smhxF(x), 



(0) 



where x — X ab I = 0. Here, r s is a scale-dependent 
interaction parameter that is renormalized as dr s /dl = 



f/4 [12j with the solution 



f s /(l + f s Z/4). Eq. 



([6]) with this ^-dependent r s is solved best by introducing 
y = ln(l + f s £/4), such that dx/dy = — 4 sinhx-F(x)- 

We start the detailed analysis of Eq. ((6]) assuming x 3> 
1. As long as x > 1 we use lim x ^oo -F'(x) = 2/37T to find 



dy 37r 



(7) 



Eqs. (O and ([7]) make two remarkable predictions: 
First, the minus signs in both equations imply that the 
Coulomb interaction reduces the barrier strength. It thus 
enhances transport through the studied scatterer, in ev- 
ident contrast with the familiar situation in ID 
Second, that enhancement is exponentially large in x 
such that for sufficiently large x the Coulomb interac- 
tion strongly impacts electron transport even at r s l <C 1. 

Both predictions have their origin in bound states. We 
find that the interaction with the Friedel oscillations due 
to propagating states decreases the conductance as in 
ID. That suppression, however, is overcompensated by 
bound states which form inside the scattering barrier at 
\x\ < a/2, but leak into the region \x\ > a/2. Electrons 
in bound states have the same probability for being to 
the left of the barrier as to the right, guaranteed by V- 
symmetry [28j. The bound states therefore straddle the 
barrier and the exchange interaction with electrons in 
those bound states opens an additional path for electrons 
to traverse the barrier. This explains the minus signs in 
Eqs. ([6]) and ([JJ). The resulting enhancement of the bar- 
rier transmission t is of order r s l and can easily exceed the 
exponentially small noninteracting t cc exp(— x), which 
implies the exponential renormalization of Xi Eq. (0- 

This exponential renormalization limits our RG anal- 
ysis to r s <C 37rexp(— x). Also, the bound state ener- 
gies e° = ±wfcj,sechx are exponentially small at x ^ 1- 
The exponential renormalization of x is thus cut-off at 
K b = kT coshx|fc'=K b / w an d it requires exponentially low 
temperatures kT <C exp(— x)v/a. Eq. ([7]) is solved by 
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(8) 



where y is evaluated at n h , y = ln[l — (f s /4) ln/t a]. The 
self-consistent solution of Eq. © with this x-dependent y 
determines the normalized conductance, which for )(>! 



takes the form G r = 8exp(— 2%)/3. In analogy with 
the Kondo problem we define the temperature scale T* 
where the transmission across the barrier becomes of or- 
der unity, G r ~ 1. At r s 37rexp(— x) we find 



kT* = - e -^/ f ° 
a 



(9) 



At wavevectors k' < k° the bound states do not renor- 
malize x anymore, but the propagating states continue to 
do so and dx/dy = 3/2 for K b > k 1 > kT/v. This now 
decreases the conductance, by a factor (l+f s x\k'= K b /4) -3 
if xlfe'=K b ^ 1- I n Fig. [5] we plot the normalized resis- 
tance GjT 1 for various parameter values. Some of them 
are within the regime f s <C 37rexp(— x), where Eq. ([5]) is 
rigorously justified, others at f s < 37rexp(— x), where Eq. 
(J5J) may still be qualitatively correct. At f s > 37rexp(— x) 
Eq. (j5J) does not have a unique solution anymore. Clearly, 
the predicted effect is large even for moderate x ~ 1 an d 
ln« b a > 1. Its observation therefore does not require 
access to large temperature intervals. 
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FIG. 2: Normalized resistance G7 1 at \ — 5 and kT < 
usechx/a, when it is renormalized by bound states, as a func- 
tion of Al = In (v sech x/kTa). 

As x decreases the scaling behavior crosses over into a 
second regime, where x 1 and F(x) may be expanded 
around y = 0. In that case we find dx/dy = — 4(10/37r — 
l)x 3 [2 9| | such that, now assuming for simplicity j(<Cl, 
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^1 + 8x2(10/3^-% 



(10) 



We conclude that the backscattering from A disappears 
asymptotically at low T, very differently from static scat- 
tering potentials v in ID. The coupling Xi however, is 
only marginally irrelevant here (dx/dl oc x 3 ,30]). It is 
renormalized more slowly than v in ID, which is relevant 
in the RG sense {dv/dl oc v) Q. The renormalization of 
r s in graphene further slows down the RG flow of x- 

Discussion: In reality the unscreened and ballistic seg- 
ments of graphene around the barrier in Fig. [1] have a fi- 
nite length L (typically limited by contacts or by the elas- 
tic mean free path) . Our predictions then require that the 



4 



entire Friedel oscillation fits into those segments, imply- 
ing an additional cut-off k l — cothx/L for the RG-flow 
and necessitating L ^> a. Due to the strong renormal- 
ization at % ^> 1, however, the observation of the effects 
shown in Fig.[5]does not require excessively long samples. 

While inelastic processes are negligible in the above 
analysis of the renormalization of x at <C 1 [l9| . 
they can significantly affect the transport between bar- 
rier and contacts [20(. There inelastic scattering can be 
neglected only if the barrier is more resistive than the 
adjacent ballistic regions, such that (fcTG r ) _1 3> r 2 L/v. 
The cut-off k l further limits the regime where our pre- 
dictions can be observed directly in the conductance to 
r s 2 G r < v/kTL < tanhx- At r 2 s G, > v/kTL the obser- 
vation of GjT 1 requires experimental differentiation from 
the inelastic background resistance, e.g. by varying \. 

The vector potential A is induced by strain in the re- 
direction if that is along the "armchair" direction of the 
graphene lattice [13] • Strain, however, also gene rates a 
scalar, "deformation" potential V < xv/a [171 ] . This 
breaks PH-symmetry, a Hartree-potential appears, and 
our predictions require V <C max{w cothx/L, kT}/r s . 
Also the bound states are strongly affected by V and 
consequently Eq. ([7]) holds only if V <C x ex P(~x) u / a i a 
tight constraint when % ^> 1. We remark that Eq. © 
(strictly valid at ^ > 1, but possibly qualitatively cor- 
rect also at x — 1) predicts strong effects like those of 
Fig. [2] down to % ~ 1, when this extra constraint dis- 
appears. To avoid a nonzero V, A can be generated 
alternatively by a pair of wires at x = 0, but different 
distances d\,d,2 — a from the graphene plane (see Fig. 
[I}. Opposite currents I and —I through such wires cre- 
ate a vector potential that falls off as 1/x 2 at a; > a and 
that has integral X = f dx A y {x) — fioI(di — c^). In our 
limit r s < 1, but I 3> 1 the deviations of this potential 
from the square- well shape Eq. (p} have negligible effects. 

Conclusions: We have shown that the Coulomb repul- 
sion between electrons can have profound and unconven- 
tional effects on transport in intrinsic graphene. General 
considerations predict singular transport across energy- 
independent ID defects without Hartree potential. In ad- 
dition we have shown that interactions in graphene can 
enhance transport through static scatterers, in surprising 
contrast with the familiar result for ID conductors. We 
have presented detailed calculations for vector potential 
scatterers that are induced by strain. There the above 
effects have exponentially large experimental signatures. 
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